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Topological insulators have an insulating bulk but a metallic surface. In the simplest case, the 
surface electronic structure of a 3D topological insulator is described by a single 2D Dirac cone. A 
single 2D Dirac fermion cannot be realized in an isolated 2D system with time-reversal symmetry, but 
rather owes its existence to the topological properties of the 3D bulk wavefunctions. The transport 
properties of such a surface state are of considerable current interest; they have some similarities 
with graphene, which also realizes Dirac fermions, but have several unique features in their response 
to magnetic fields. In this review we give an overview of some of the main quantum transport 
properties of topological insulator surfaces. We focus on the efforts to use quantum interference 
phenomena, such as weak anti-localization and the Aharonov-Bohm effect, to verify in a transport 
experiment the Dirac nature of the surface state and its defining properties. In addition to explaining 
the basic ideas and predictions of the theory, we provide a survey of recent experimental work. 
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I. INTRODUCTION 

Topological insulators (TI) are bulk insulators with 
protected metallic surface states as a result of the topo- 
logical properties of the bulk electronic wavefunctionsP^ 
For a three-dimensional topological insulatorp' this 
metallic state is a two-dimensional electron gas with 
many special features such as spin-momentum locking 
and a robustness to localization by disorder. While sev- 
eral of these features have been observed with surface sen- 
sitive probes such as angle-resolved photoemission spec- 
troscopy, much of current experimental focus is aimed at 
demonstrating these and other surface state properties in 
transport measurements. 

The main goal of this review article is to explain 
how the Aharonov-Bohm and other magnetotransport 
effects are manifested in topological insulator surface 
states and to summarize recent experimental and theoret- 
ical progress towards their observation. The Aharonov- 
Bohm effect can be utilized as a fundamental probe of 
how the quantum phase of an electronic wavefunction is 
sensitive to magnetic flux through the gauge invariance 
of the Schrodinger equation coupled to electromagnetic 
fieldsP^It is possible to give self-contained explanations 
of how the features of the surface state, which is mod- 
eled by a single massless Dirac cone in the simplest case, 
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lead to unusual (as compared to in a traditional two- 
dimensional electron gas) magnetotransport behavior in 
a variety of experimentally relevant situations. This ap- 
proach leaves out only the connection between bulk wave- 
functions and surface electronic states, which requires a 
little bit of mathematical background and has been ex- 
plained several times, for example in the reviews cited 
above PHI 

In the remainder of the introduction we discuss the 
basic properties of topological insulator surface states, 
concentrating on what makes them different from the 
two-dimensional electron gas (2DEG) in either conven- 
tional semiconductor heterojunctionsP or graphenePESl 
The main properties of the most common topological 
insulator materials are summarized and the relevant el- 
ements of q uan tum transport including the Aharonov- 
Bohm effectP^ and localization theorjffiffi^ is given. 
While this review is mainly focused on the transport 
properties of 3D topological insulators, we comment at 
various points on the two-dimensional topological insu- 
lator or quantum spin Hall state, 14 as observed first in 
(Hg,Cd)Te quantum wellsf^ This phase is is interesting 
in itself and also illuminates some aspects of the three- 
dimensional behavior. 

The main text consists of four sections. In section ITT1 we 
discuss the longitudinal conductivity which in the ideal 
case of insulating bulk should show ambipolar Hall effect 
and minimal conductivity when the chemical potential is 
tuned through the Dirac point. The main theoretical re- 
sult covered is the absence of localization and the accom- 
panying flow to the symplectic metal phase. Section |III| 
addresses magnetic field induced quantum oscillations 
such as the Shubnikov-de Haas (SdH) oscillations. We 
focus in particular on the signatures of the Berry phase 
of the Dirac fermion in the SdH signal. The quantum Hall 
effect is briefly discussed. Magnetic flux effects on quan- 
tum transport, such as weak anti-localization (WAL) and 
Aharonov-Bohm (AB) oscillations, are the subject of sec- 
tion |IV| WAL is a useful probe of 2D transport but is 
insensitive to the Berry phase. In contrast, the AB oscil- 
lations can in principle be used to infer the presence of 
a nontrivial Berry phase. Finally, in section [V] we collect 
some related problems and future directions. We end the 
review with a summary and conclusion. 



A. Surface states of three-dimensional topological 
insulators 



In this section we introduce some key definitions and 
explain how the surface state of a topological insulator 
differs from a conventional two-dimensional electron sys- 
tem (cf. figure [TJ. Any strictly two-dimensional metal 
with time-reversal symmetry has a Fermi surface that 
consists of an even numb er of "sheets" (closed curves) 
once spin is included! 16 * 17 ! In the simple case of no spin- 
orbit coupling, the two sheets are degenerate; they move 
apart when spin-orbit coupling is included, but at every 



value of the Fermi energy the Fermi surface still consists 
of an even number of closed curves. The surface state 
of a topological insulator is not strictly two-dimensional 
in the sense that it consists of a boundary between two 
different three-dimensional bulk states, one of which is 
frequently the vacuum. It has an odd number of Fermi 
surface sheets, and in the simplest case that odd number 
is 1. 

As an example, consider the model linear dispersion 
relation obtained with the Hamiltonian 



H = V(p X (Jy -Py(J X ), 
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where v is the Fermi velocity, p the momentum, and a are 
the Pauli matrices. Another commonly used Hamiltonian 
H = v p • cr differs only in the interpretation of the sigma 
matrices and their relation to the real spin. For experi- 
mentally relevant materials the details can be more com- 
plicated: the velocity is not simply scalar but depends on 
the direction GUMD the spin p olarization may be reduced 
from its maximal value) 21 " I there may be more sheets of 
the Fermi surface^ (i.e., an odd number larger than 1) 
and nonline arity of the spectrum shows up away from the 
Dirac point pOUl Details of the surface geometry can also 
influence the electronic and spin structure I2SH2QI However, 
the Aharonov-Bohm effect and other magnetotransport 
phenomena discussed below are generally independent of 
these details. The "Dirac cone" energy-momentum rela- 
tion described by Q is shown in figure [T 

Time-reversal symmetry implies that the state at mo- 
mentum p must have spin direction opposite to that at 
—p. The spin direction therefore precesses as the elec- 
tron momentum moves around the Fermi surface. This is 
often referred to as spin-momentum locking. This case is 
different from a "half-metal" where the spin polarization 
is constant and time-reversal symmetry is broken. As an- 
other comparison, adding the quadratic term p 2 /2m to 
the above Hamiltonian gives a Hamiltonian commonly 
used to describe quantum wells with Rashba spin-orbit 
coupling!^! At every energy, there are either zero or two 
sheets of the Fermi surface (see figure [T| . 

In 3D, the bulk wavefunctions of a perfect crystal are 
characterized by four topological invariants that take val- 
ues in Z 2 , three "weak" and one "strong"|22lEH] The latter 
is the one of greatest interest both experimentally and 
theoretically, and we will use the term "topological insu- 
lator" to refer to strong topological insulators where this 
invariant is nonzero. The existence of an odd number of 
Fermi surface sheets is a consequence of the odd value of 
the strong topological invariant. There are also some in- 
teresting features if the strong index is zero but one of the 
weak invariants is nonzero. Such "weak topological insu- 
lators" ( WTI) have an even number of Dirac cones on the 
surface and can be viewed as layered versions of the two- 
dimensional topological insulator or quantum spin Hall 
state. The surface states of WTI's are in principle less 
stable to either bulk or surface disorder since there ex- 
ists perturbations that gap out the surface. It turns out, 
however, that when this perturbation is zero on average 
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FIG. 1. Comparison between the Dirac dispersion (left top 
and bottom) and the dispersion of a 2DEG with Rashba spin- 
orbit coupling (right top and bottom). Both dispersions are 
rotationally symmetric around the energy axis, as shown by 
the plot of the Fermi surface in the lower panel. The Fermi 
surface consists of an odd and even number of closed curves 
(sheets) respectively. The filled and crossed circles and ar- 
rows denote the spin direction of the corresponding eigen- 
state. Due to time-reversal symmetry, scattering from a spin 
up state at k to a spin down state at — k is forbidden (crossed 
out dashed arrowed lines). Backscattering is therefore com- 
pletely absent for the Dirac dispersion. In the Rashba case, 
scattering between the branches with the same spin state is 
allowed (dashed line without a cross), and backscattering can 
take place. 



as in a random environment, the surface state is robust 
against localization. We discuss this in more detail in 
section IV Al 

The low energy electronic structure of graphene is also 
described by two Dirac cones (ignoring the spin which 
does not play an important role in most graphene ex- 
periments). Graphene, however, has a time-reversal T 
with T 2 = +1 and is therefore in a different symmetry 
class from TPs which time-reversal symmetry satisfies 
T 2 = — 1. In the absence of intervalley coupling an ef- 
fective time-reversal symmetry with T 2 = — 1 emerges 
in graphenej^ and the physics is that of a single Dirac 
cone. It is useful to keep this in mind since several re- 
sults originally obtained for graphene are relevant to TI 
transport. 

The two-dimensional topological insulator or quantum 
spin Hall state has the same locking of spin and momen- 
tum in its edge state as in the surface state described 
in 0: electrons moving one way along the edge have a 
certain spin orientation which is opposite that of the spin 
of the electrons moving in the reverse direction. There is 



a single branch of edge excitations moving in each direc- 
tion, unlike in an ordinary quantum wire, which has two 
branches (spin-up and spin-down). 

One simple difference in surface state transport be- 
tween the 2D and 3D topological insulators can now be 
explained, and this will also help convey the importance 
of time-reversal symmetry. Time-reversal symmetry im- 
plies that every spin-half cigenstate is degenerate in en- 
ergy with, and distinct from, its time-reversal conjugate 
(the state obtained by reversing the direction of time) . As 
a result, every energy eigenvalue in a time-reversal invari- 
ant system of independent electrons is at least two-fold 
degenerate; these degenerate pairs are called Kramers 
pairs. Integer-spin particles can be equivalent to their 
time-reversal conjugates and there need not be such de- 
generacies. Now consider perturbing the original Hamil- 
tonian. The robustness of Kramers pairs necessitates 
that any time-reversal invariant perturbation, such as po- 
tential scattering, has a zero matrix element between the 
two states of a pair, as otherwise it would split the pair. 

As a consequence of the robustness of the Kramers 
pairs, elastic scattering at the edge of a 2D topological 
insulator disappears at low energy (in the gap), because 
the two available states belong to the same Kramers pair. 
The corresponding scattering from spin-up to spin-down 
is still forbidden in the case of an ordinary wire, but the 
scattering process that does not flip the spin is allowed 
and eventually leads to localization by disorder. In gen- 
eral, an even number of Kramers pairs of edge modes 
will localize, while an odd number will lose pairs until 
a single pair is left which cannot be localized. At low 
voltage and temperature, transport is effectively ballistic 
because backscattering disappears, although the correc- 
tions to the quantized conductance e 2 /h are expected to 
be power-law rather than exponential as in the quantum 
Hall caseP 

The same Kramers protection exists at the surface of a 
3D topological insulator but is less powerful as now there 
are allowed scattering processes that do not violate the 
Kramers theorem. Scattering at any angle other than 180 
degrees is allowed, and indeed Fourier transforms of STM 
measurements^ show the vanishing amplitude of perfect 
backscattering. Because there are still allowed scatter- 
ing processes at leading order, unlike in the 2D case, the 
low-temperature fate of conduction at the surface of 3D 
topological insulator requires more thought and is dis- 
cussed below in the context of weak localization theory. 



B. Properties of topological insulator materials 

In this review, we will focus on the basic physics of the 
topological insulator phase revealed through Aharonov- 
Bohm measurements and other magnetic effects on trans- 
port. However, in order to understand experiments, 
it seems useful to provide a few notes on the mate- 
rials studied in current experiments, even though ma- 
terial improvements are rapid. In two dimensions the 
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first demonstration of the theoretically predicted "helical" 
edge state was in (Hg,Cd)Te quantum wells Recently, 
experiments were reported showing evidence for helical 
edge channels in InAs/GaSb quantum wellsP^ There are 
theoretical proposals to realize the phase in other sys- 
tems, e.g., when heavy atoms are adsorb ed on graphene 
to increase the spin-orbit couplin g* 39 ! 40 ! and in strained 
graphene in the presence of interactions.^] 

In the remainder of this section, we concentrate on 
the 3D state, where there are more materials and exper- 
iments. Bi-Sb alloys were the first materials studied for 
topological insulator behavior but the high level of al- 
loy disorder and the complicated surface Fermi surface 
(with 5 band crossings along the cut studied) have led 
to their being superseded by other materials, in partic- 
ular the semicon ductor s Bi2Se3, Bi2Te3, Bi2Te2Se, and 
variations thereofP^EH 

Bi2Se3 is the 3D TI material that has been most inves- 
tigated experimentally.^ It has a trigonal unit cell with 
five atoms, and can be pictured as having five layers Se- 
Bi-Se-Bi-Se; the central Se layer is clearly inequivalent 
to the outer two layers, and there is a good cleave plane 
between the van der Waals-bonded first and last Se lay- 
ers. The bulk bandgap is approximately 0.3 eV. The sur- 
face Dirac velocity depends somewhat on the energy and 
the direction (there is a significant hexagonal distortion 
except very close to the Dirac point ps3) but an approxi- 
mate value useful for theoretical estimates is 4 x 10 5 m/s. 
These numbers are consistent with estimates from GW- 
improved DFT calculations! 4 - Estimates of the surface 
state spin polarization from photoemission and numeri- 
cal calculations ra nge fro m 60% of the spin half maximum 
up to nearly 100%J2U22I 

Bi2Te3, also shown to be a topological insulator 
has been studied for many years as a practical room- 
temperature thermoelectric material. It has the same 
structure as Bi2Se3 and a smaller bulk bandgap of 0.15 
eV. Its thermoelectric utility can be understood from 
the rule of thumb that the ideal operating temperature 
of a thermoelectric semiconductor is about one-fifth of 
the bandgap fcsT w E g /h. A significant power fac- 
tor (product of electrical conductivity, temperature, and 
thermopower squared) depends on having an appreciable 
number of thermally excited carriers. Bi2Te2Se ("BTS") 
has been studied very actively because crystals can be 
grown with bulk conductivity order s of magnitude lower 
than either Bi2Se3 or Bi2Te3.^E^ The Se atoms go in 
the middle layer of the five-layer structure. As opposed 
to Bi2Se3, which has the Dirac point in the bulk gap, the 
Dirac point in Bi 2 Te3 is buried deep in the valence band. 

Considerable effort is going into finding 3D topolog- 
ical insulators and related phases in other materials 
families, whether to decrease the bulk conductivity or 
to find ordered phases that combine topological order 
with another type of order (e.g., antiferromagnetisnPH or 
superconductivity.^) The magnetotransport effects de- 
scribed below may be useful in identifying materials in 
the topological insulator phase when angle-resolved pho- 



toemission or other methods are impractical. In addi- 
tion to the Bi based materials, straining 3D HgTe which 
is nominally a semimetal, opens up a gap and real- 
izes a TlP^ /3-Ag2Te has also been used in transport 
studies JS2152I 

In early transport experiments the conductance was 
dominated by the bulk. The cleanest signatures of quan- 
tum interference from the surface state have been ob- 
tained in thin films and nanowires. These are cither 
epitaxially grown or obtained with mechanical exfolia- 
tion. Due to the easy cleave plane, the thickness is gen- 
erally a multiple of quintuple layers with each quintu- 
ple layer about 1 nm thick. In the ultrathin limit the 
tunnel coupling between the top and bottom surface is 
sufficiently large to open up a sizable gap and make the 
film insulating. 1 " 56 l This gap is obse rved experimentally 
in films a few nanometers thickP^^ and its rapid decay 
to zero with increasing thickness gives a direct measure- 
ment of the surface state's penetration into the bulk that 
can be compared with photoemission results.^ Insulat- 
ing behavior in transport has also been observedP^HISl 

With increased surface mobility the possibility of re- 
alizing correlation physics opens up. One direction is 
fractional quantum Hall physics at the surface; indeed 
some features are observed in transport data that are con- 
jectured as possibly indicating incipient fractional Hall 
states at the surface!®! A thin film of topological insula- 
tor in a magnetic field is analogous to a quantum Hall 
bilayerP^l These bilayers have been a fertile system for 
studies of many-body physics, as the proximity of the 
two layers enables interlayer correlated phases with re- 
markable properties when the Coulomb interaction be- 
tween the layers is significant. However, these correla- 
tion effects may be hidden by the hybridization gap in 
topological insulator thin films. 



C. Diffusive transport and the symplectic metal 

In the Altland-Zirnbauer symmetry classification,^ 
the topological insulators we are interested in here be- 
long to the symplectic class (All in the Cartan notation) . 
This symmetry class is characterized by the presence of 
a time-reversal symmetry. The time-reversal operator T 
commutes with the Hamiltonian and satisfies T 2 = — 1. 
A two dimensional electron gas (2DEG) with strong spin- 
orbit coupling is also in the symplectic class, as well as 
graphene in the absence of intervalley scattering J^l When 
the electronic transport is diffusive many of its charac- 
teristics are independent of the topological properties of 
the underlying system. We will refer to this state as the 
symplectic metal. The symplectic metal is characterized 
by weak anti-localization (WAL). Because of WAL the 
symplectic metal is, in the renormalization group sense, 
a stable fixed point. In this section we review the var- 
ious characteristic quantum interference phenomena of 
the symplectic metal. We limit the discussion to two ter- 
minal transport. This ideal setup can for example model 
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FIG. 2. A schematic of a typical two terminal transport setup 
(left), and the ideal modeling with a set of incoming and out- 
going modes (right). A current is injected through one lead 
and extracted through the other. In the ideal model, the 
current is assumed to go only through the surface and not 
through the bulk. We allow for the possibility that the modes 
are different in the left (L) and right (R) lead, though in prac- 
tice the modes in the two leads are often related to each other. 
For example, in the Dirac case |3| \n) l = |7*(— up to a 
possible phase. 



transport in one side of a large TI surface or conductance 
through a TI nanowire as schematically shown in figure [2] 

The presence of a time-reversal symmetry crucially af- 
fects many of the transport properties of the symplec- 
tic metal. Before discussing the perhaps more intuitive 
path picture of the interference phenomena, we consider 
how these directly result from the symmetry constraints 
imposed by time-reversal symmetry on the scattering 
matrix describing the two terminal transport. In par- 
ticular, we obtain absence of backscattering, 6 ^ discuss 
under what condition a perfectly transmitted mode is 
realized,^ and give the connection between the time- 
reversal symmetry and weak anti-localization. 



1. Quantum transport and time-reversal symmetry 

The scattering matrix describing the two terminal 
transport in the presence of time-reversal symmetry with 
T 2 = — 1 can be chosen to be antisymmetric. Before dis- 
cussing the consequences of this antisymmetry we give a 
short derivation of this fact, following Ref. [68] 

Consider a two terminal setup as in figure [2] The left 
and right leads are metallic contacts which host a large 
number of incoming modes, denoted by \ti)l and \n) R 
respectively. In a topological insulator these modes are 
the properly normalized eigenstates of the Dirac Hamil- 
tonian . For example, if we assume periodic boundary 
conditions in the transverse y direction, the modes can 
be written 



Hi 



\/2 \ i 



a ikx+iq n y 



(2) 
(3) 



Here q n = 2-KnjW are the discrete transverse momenta 
with W the transverse width of the sample, and the Fermi 
energy Ep — hvy/k 2 + q 2 . We have assumed k 3> q n 



such that we can ignore the momentum dependence of 
the spinors. This is the relevant limit in the leads, which 
being metallic are highly doped These modes carry 
unit current (n\ <j x \ti)l = — {n\ u x |n)fl = 1. The outgo- 
ing modes are the time-reverse of the incoming modes 
\7~n) l and \Tn) R , since time-reversal reverses the direc- 
tion of motion. Importantly, \l~n) has the opposite spin 
to \n). 

A solution to the Hamiltonian describing the sys- 
tem, subject to the boundary conditions imposed by the 
metallic leads at x = and x — L, can be written 

Y. n c n,L\n) L + d n , L \Tri) L , x<0, 
Z)„ Cn,R\n) R + d n . R \Tn) Rl x > L, ! 4) 
>Af), 0<x<L. 

Since is a solution to the Schrodinger equation, the 
incoming coefficients cl and c R are linearly related to the 
outcoming coefficients dj, and d R through the scattering 
matrix S 



d L 
d R 



(5) 



Assuming an equal number of modes N on the left and 
the right, the x 2N scattering matrix has the block 
diagonal form 



S = 



r t' 
t r' 



(6) 



where r nm is the probability amplitude of reflection of 
mode \m)i l into |TVi)l, and the other block matrices 
have a similar interpretation. For our current purpose 
we do not need to know the exact form of the wavefunc- 
tion \iPm) i n the sample. In an actual calculation the 
scattering matrix is obtained from a knowledge of \i/jm)- 
The scattering matrix gives the two terminal conductance 
through the Landauer formula 



e 2 e 2 
G = —tith = — tr (1 - rV), 
h h 



(7) 



where the second equation follows from current conser- 
vation S^S = 1. For ease of notation we also introduce 
the unitless conductance g — G/(e 2 /h). The conduc- 
tivity of a sample of width W and length L is given by 
a = GL/W. 

Because of time-reversal symmetry, if \ip) is an eigen- 
state of the Hamiltonian so is \T^p). This solution, fur- 
thermore, is orthogonal to \ip) since 



(7>| t/j) = (T 2 ^| Ti/j}* = - (TV>| = o. 



(8) 



In the first equality we used the antiunitarity of T, then 
that T 2 = — 1. \Ttfj) is therefore an independent solution 
and all energy eigenvalues are doubly degenerate: the 
celebrated Kramers degeneracy. As a consequence 



-n,L 



\Tn) L - d* L \n) L , x < 0, 



|TV> = <T, n <,R\Tn) R -dt R \n) R , x > L, (9) 



\Ti>u), 



< x < L, 
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is also an allowed scattering state. By inspecting Q we 
obtain 



" 1 



(10) 



Comparing with Eq. (|5| one concludes that time-reversal 
symmetry requires 



-S. 



(11) 



with T denoting the transpose. The antisymmetry of the 
scattering matrix gives rise to the absence of backscat- 
tering, the presence of a perfectly transmitted mode, and 
weak anti-localization, as we will now explain. 

The absence of backscattering is simply the fact that 
the diagonal elements r nn — r' nn = 0. This statement 
holds true both in a TI and a 2DEG, since it only requires 
the presence of a time-reversal symmetry. However, since 
r nn is reflection of mode \n) back into \Tn), which has the 
opposite spin, backscattering that does not flip the spin 
is allowed in the 2DEG (cf . figure [T]) . There is no such 
state in a TI and therefore backscattering is completely 
absent. 

Due to the antisymmetry of r, the eigenvalues of r'r, 
and by unitarity also the transmission eigenvalues of tH, 
come in degenerate pairs. This is the Kramers degen- 
eracy of transmission eigenvalues. Furthermore, if the 
number of modes N is odd, there is necessarily at least 
one transmission eigenvalue that is equal to unity, since 
det(r) = det(-r T ) = (-1)"^ det(r) = -detr = 0. This 
is the perfectly transmitted mode discussed by Ando and 
SuzuuraJ^and it will play a central role in our discussion 
below. 

An odd number of modes can never strictly be re- 
alized in an inherently 2D system with time-reversal 
symmetry! 16 * 17 ! The absence of certain couplings can, 
however, reduce the system to an effective one with an 
odd number of modes. An example of this is graphene 
in the absence of intervalley coupling. 35 In contrast, a 
TI can intrinsically have an odd number of modes and 
therefore host a perfectly transmitted mode. 

Weak anti-localization is the first order in 1 / g quantum 
correction to the classical Drude conductance. To obtain 
WAL from the scattering matrix we assume the elements 
of S to be randomly distributed Gaussian variables 70 
with 



1 finm 

2N- 1 



(12) 



where the angular brackets denote average over the dis- 
tribution of the elements of the scattering matrix. The 
delta function is needed to satisfy the antisymmetry con- 
dition (111 and the denominator is obtained from the 
unitarity of S. Using this in the Landauer formula ^ 



g = N- 



N 2 -N 
2N - 1 



N 
~2 



O 



N 



The first term is the classical conductance. It takes the 
value N /2 since each mode is equally likely to be trans- 
mitted as being reflected, due to the randomness of the 
scattering matrix. The second term is a positive enhance- 
ment of the conductance due to quantum interference, 
namely weak anti-localization. This term is absent in 
the absence of time-reversal symmetry, as is readily ver- 
ified. A time-reversal breaking pertubation, such as a 
magnetic field, will therefore decrease the conductance. 
The fact that the first quantum correction is positive re- 
flects the stability of the symplectic metal phase to weak 
disorder. This correction is perturbative (here in 1/N) 
and independent of topology and is therefore the same 
for TI surfaces and a 2DEG. The case of strong disorder 
and small conductance is discussed in section III Bl 

This argument is instructive in that it shows the re- 
lation between WAL and time-reversal symmetry. The 
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assumption ( 12 1 is however only strictly valid when the 
two leads are connected by a quantum dotP^ We are 
interested in the case of a two dimensional sample con- 
necting the two leads. This will be discussed in the next 
section. 



2. Weak anti-localization and Berry phase 

When a Dirac fermion traverses a loop in space, the 
spin rotates by 2ir due to the spin-momentum locking. 
The wave function, being a spinor, acquires a phase of it 
which can alternatively be considered as a Berry phase 
induced by the Dirac point. This phase affects quantum 
interference and in particular changes the constructive 
interference of spinless fermions that gives rise to weak 
localization into destructive interference and WAL. 

Studies of weak (anti)-localization date back a couple 
of decades and th e phys ics is by now well understood. A 
number of reviewaHEEl an d textbooks (our discussion is 
of similar flavor as^) exist that discuss the basic phe- 
nomena. The importance of spin-orbit coupling and the 
resulting WAL correction was first derived by Hikami, 
Larkin and Nagaoka.E21 An intuitive picture has been 
given by Bergmann.^ The details of the derivation for 
Dirac fermions are slightly different, but the final an- 
swer for the WAL correction is the same, as shown by 
Suzuura and AndcES and McCann et alP^ in the context 
of graphene. This is to be expected since the two systems 
are in the same universality class. Nevertheless, we in- 
clude here a short introduction to the topic, focusing on 
the broad physical picture and avoiding detailed formal- 
ism. This should serve as a guide to the literature and 
to make this review more self-contained. In particular, 
we will need some of the results discussed in this section 
in our later survey of TI transport experiments. In the 
context of TI's the WAL correction has been discussed 
further theoretically by several authors.^^ 

As in the last section we focus on the two terminal 
setup and the conductance in terms of the scattering 
matrix. Reflection and transmission amplitudes can be 



FIG. 3. Classical paths that interfere to give rise to weak 
anti-localization, a) The path a and its time-reverse a, which 
return to the time-reversed mode in the lead, interfere de- 
structively leading to the absence of backscattering. Paths f3 
going from one lead to the other do not interfere with their 
time-reversed partner j3 since these paths start in the oppo- 
site lead. However, b) paths 7 (7) that have a single (avoided) 
crossing do interfere and including these paths is crucial to ob- 
tain a current conserving theory. In addition to the 7 paths 
shown, there will be similar paths that return back to the 
same lead after the crossing. 



written as a sum over all possible paths connecting the 
two leads (see figure [3]) . 

r„ = E 4e« s »/ s , (14) 

turn = EpAf,e iSf,/h - (15) 

Here a denotes paths that enter the sample from mode 
|m) l in the left lead and exit through mode \Tn) l in the 
same lead. Similarly, /3 denotes paths that start in mode 
\m) l in the left lead and exit through mode \Tn)n on the 
right. S a is the classical action of path a and A a con- 
tains the stability amplitude, no rmaliz ation and possible 
geometric phases of the path a. 82 83 In particular, any 
phase acquired from rotation of the spin degree of free- 
dom will enter through A a . Generally, the effect of the 
spin rotation is captured by a matrix valued amplitude^ 
but this is not needed for our current purpose due to the 
spin-momentum locking. 

In terms of reflection amplitudes, the conductance ^ 
is given by 
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N ■ 



i(S a -S a ,)/h 



(16) 



When averaging over disorder, the quickly oscillating 
exponential term will generally average to zero, unless 
S a ~ S a '. The classical Drude conductance is obtained 
by including only the diagonal terms a' = a. In addition, 
interference terms where oi = a with a the time-reverse 
of a survive the disorder average since by time-reversal 
symmetry S a = S a - The amplitudes A a and A a have 
the same absolute value, but have a phase difference of 
7r, and thus A a = — A a . This is due to the Berry phase 
picked up by the 2ir relative rotation of the spin between 
the two paths, see figure [4] Therefore the total contribu- 
tion of these two paths to the sum in ( 16 1 is 



FIG. 4. The spin rotation along the paths a and a. Since the 
paths enter with opposite sign in the interference term the 
total spin rotation is 2n, resulting in the important minus sign 
that causes destructive interference in the reflection amplitude 
and the absence of backscattering. 



This is the absence of backscattering r nn — obtained 
from symmetry in the last section. Compared to the clas- 
sical value 2|^4 Q | 2 , the probability to reflect back into the 
same lead is reduced and the conductance is enhanced. 
Note that this destructive interference only affects the di- 
agonal elements of the matrix r of reflection amplitudes. 

This is not the full story, as one infers from the fact 
that such paths do not affect the conductance when writ- 
ten in terms of transmission amplitudes 
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Y,Af,A}, t 

0,0' 



i(S e -Sp,)/H 



(18) 



Indeed, if j3 goes from the left lead to the right lead, the 
time-reversed path j3 goes from the right lead to the left 
lead and therefore does not enter the sum (18). With no 



\A a \ 2 + \A a \ 2 + 2^(A a A* a )=0. 



(17) 



change to the transmission amplitudes as compared with 
the classical value, current is no longer conserved and the 
scattering matrix is not unitary if we include only these 
paths. To get a consistent, current conserving theory, we 
need to include additional paths. These are the paths 7 
that have a single self crossing (figure |3| , and the corre- 
sponding path 7 which is almost the same except for an 
avoided crossing at the crossing point. 85 The two paths 
therefore traverse the loop in the path in opposite di- 
rection. The sign of the correction to the amplitude is 
determined by the phase picked up in the loop and the 
details of the encounter at the crossing.^! Instead of go- 
ing through the derivation, we can simply infer the sign 
of the correction from the absence of backscattering and 
the need to obtain a unitary scattering matrix. In addi- 
tion to the paths 7 that end up in the opposite lead after 
leaving the loop, there are paths that instead return to 
the same lead and therefore affect reflection probabilities. 
The interference of the loop is independent of which lead 
the particle ends up in after leaving the loop. Thus, the 
correction from these loops has the same sign for reflec- 
tion and transmission probabilities. This correction is 
therefore necessarily positive to compensate the negative 
contribution due to the absence of backscattering. Trans- 
mission is enhanced by the quantum interference giving 
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FIG. 5. The different parts of the interference between the 
paths with a single (avoided) crossing in figure [3] can be in- 
terpreted as objects that are obtained in diagrammatic calcu- 
lations. The initial and last part where the paths propagate 
together corresponds to the diffusion (top) while the propa- 
gation in the opposite direction in the loop corresponds to 
the cooperon (middle). The crossing where the paths switch 
between propagating in the same direction and opposite di- 
rection corresponds to the Hikami box in the diagrammatic 
language (bottom). Here we have shown the dressed Hikami 
box which is obtained since all the terms on the right are 
of the same order. The crossing on the left is topologically 
equivalent to crossing as plotted in figure [3] 



rise to WAL. In this sense absence of backscattering and 
WAL are integrally related through current conservation. 

In the next step we focus on the quantum corrections 
to the conductance as interference between the paths 7 
and 7 as depicted in the diagram of figure [3]d. Instead of 
having the diagram denote only two paths we reinterpret 
it as a sum of many paths. Decomposing the diagram into 
parts, each part has an interpretation which corresponds 
to an object that is commonly obtained in diagrammatic 
calculations of the conductivity. This connection between 
the path picture of WAL and the diagrammatic calcula- 
tion aids the intuitive understanding of the latter. 

To that end, we decompose the diagram of figure (SJd 
into three main parts. The first part represents the sum 
of all possible pair of paths that travel together from 
the lead towards the crossing point. This includes paths 
that are scattered by a single impurity, two impurities, 
and so on. Since they travel together, the two paths 
scatter of the same set of impurities in the same order. 
This is equivalent to the "diffusion" in the diagrammatic 
language and naturally satisfies a diffusion equation. Es- 
sentially, a density disturbance is injected at the lead and 
diffuses towards the crossing point. The second part is 
the loop, which the two paths traverse in opposite direc- 
tion. The paths again scatter of the same set of impuri- 
ties, but now in the opposite order. In the diagrammatic 
language this corresponds to the "cooperon". When time- 
reversal symmetry is preserved, it also satisfies a diffusion 
equation. The last part is the crossing point, where the 
paths switch from traveling together to traveling in the 
opposite direction. This corresponds to the Hikami box 
in the diagrammatic language. 

To find the correction to the conductivity we need to 
estimate the number of paths with a single self-crossing. 
This is equivalent to calculating the probability for the 



cooperon to return back to the crossing point. Since it 
satisfies a diffusion equation, in time t this probability 
is proportional to l/t. It does not matter at what time 
the crossing happens, as long as the time is smaller than 
the phase coherence time after which the paths are no 
longer phase coherent. Integration over time gives 
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(19) 



The lower bound is given by the mean free time r, below 
which the motion is ballistic, and the dephasing length 
l\ = Dt^ with D the diffusion constant. In the last 
equality we have assumed the temperature dependence of 
the dephasing to be given by a power law ~ T~ p with 
p > 0. Tt is defined by the relation t /t = (T/T t )- p . 
We have not explicitly kept track of the prefactors, since 
our main goal here is to give a phcnomcnological under- 
standing of the origin of the logarithmic dependence of 
the WAL correction. Details can be found for example 
in Ref. ED 

We have described a current conserving theory of WAL 
that takes into account paths with a single self crossing. 
In principle, there are additional paths with a larger num- 
ber of crossings and these also give arise to interference 
correction. One can show that including only the single 
crossing paths is equivalent to assuming kp£ ^> 1, with 
kp the Fermi momentum and £ the mean free path.^ 
This is the same condition required for the validity of 
most diagrammatic calculations. More crossings give 
higher order corrections in l/kp£, and are eventually re- 
sponsible for localization (see e.g. Ref. [571 and references 
therein for a description of such higher order paths in 
chaotic transport). 

In theoretical work on quantum transport of Dirac 



fcrmions, the WAL correction (19) is often used to iden- 
tify the symplectic metal. For example, in section |II B| 
on the absence of localization and section IV Al on trans- 
port in WTI's, a logarithmic dependence of the conduc- 
tance on the system size (replacing the phase coherent 
length £tp in finite systems) with a slope e 2 /irh is ob- 
served. Experimentally, is is often easier to observe WAL 
by applying a magnetic field. The magnetic field breaks 
time-reversal symmetry and changes the interference of 
the cooperon such that eventually the WAL correction 
disappears. This is the subject of the next section. 



3. Aharonov-Bohm and Altshuler-Aronov-Spivak 
magnetoconductance oscillations 

In the presence of a magnetic field perpendicular to the 
2D motion, the paths a, j3, 7 pick up the Aharonov-Bohm 
(AB) phased 



»f L A - dr 



(20) 
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with A a vector potential of the magnetic field B = V x 
A, and the integration is along the path a. If the path is a 
closed loop, the integral is equal to the flux cf> through the 
loop and the Aharonov-Bohm phase becomes 2-K<f>/(f)Q, 
with cf>o — h/e the magnetic flux quantum. The AB phase 
can both lead to periodic oscillation of the conductance 
and destruction of WAL. 

We assume the magnetic field is weak enough not to 
affect the classical motion on the time scale of the mean 
free time r, i.e. that lo c t <C 1 with oj c the cyclotron fre- 
quency. The diffusion constant, and consequently the 
Drude conductivity, therefore remains unchanged. In the 
diffusion the two paths pick up the same phase which is 
then canceled. In the loop (cooperon) the two paths pick 
up the opposite phase. This modifies the diffusion of the 
cooperon which now satisfies a diffusion equation with 
p > p — 2eA. Solving the modified diffusion equation, 
one finds that the return probability is proportional to 
B/<j) sinh(AnDtB/(j) )^The WAL correction takes the 
form 
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where tb = h/(4eDB) and VP is the digamma function. 
The characteristic magnetic field strength needed to de- 
stroy the phase coherence of the cooperon, and thereby 
the WAL, corresponds to a flux quantum through the 
loop. Since Tg > t the first digamma function is often 
replaced by its asymptotic form, resulting in the Hikami- 
Larkin-Nagaoka expression^ 
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We have introduced the factor a which is simply equal 
to 1/2 here, but will be useful in comparing with experi- 
ments below. 

A magnetic field can also give rise to periodic oscil- 
lations in the conductance when the surface is multiply 
connected. This is realized for example in TI nanowires, 
where if the bulk is insulating the surface acts as a hollow 
metallic cylinder. Applying flux <j) along the wire intro- 
duces a phase 2imij) / '</>o to each path that loops around 
the cylinder n times. This modifies the WAL correction 
in a sample of width W and length L to be 
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up to an unimportant constant. The magnetoconduc- 
tance oscillates with a period of </>o/2. The factor of two 
originates in the interference between clockwise and anti- 
clockwise circulating paths that each pick up a flux </>. 
This type of oscillations where originally discussed by 
Altshuler, Aronov and Spivak in the context of metal- 
lic cylinders,^ and experimentally observed by Sharvin 
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FIG. 6. In the field theory of diffusion, the non-linear sigma 
model, the propagators of the fields correspond to the diffu- 
sion and cooperon (top, cf. also figures [3] and [5| . The four 
field interaction vertex coincides with the Hikami box (bot- 
tom). Higher order interaction terms are also present (not 
shown). 



and Sharvin!^! For this reason they are often referred 
to as AAS oscillations in the literature (for a review see 
Ref. [6]). Away from the symplectic metal limit the mag- 
netoconductance of the TI nanowire can realize robust 
oscillations with period </>o, double that of the AAS os- 
cillations. This is discussed in section HVBI 



4- Universal conductance fluctuations 

Weak anti-localization is a quantum correction to the 
average conductance. For a given sample, as external pa- 
rameters are varied, the conductance fluctuates around 
the average due to quantum interference. The ampli- 
tude of the fluctuations is universal (independent of mi- 
croscopic parameters) and of the order of e /h. These 
are the universal conductance fluctuations (UCF). Like 
WAL, the UCF amplitude is determined by the symme- 
try class and is a characteristic of the symplectic metal 
that is independent of topology. The UCF can be under- 
stood in terms of the path picture of the last section. We 
will not discuss this here and instead refer for example 
to Ref. I7T1 for details and further references. 

In the context of Dirac fermions the UCF was stud- 
ied theoretically in Refs. UCF has been observed 
experimentally in TI's in several experiments and its tem- 
perature dependence used to estimate the phase coher- 
ence length [93H5T1 Anomalously large conductance fluc- 
tuations were observed in large Bi2Se3 crystals,^ pre- 
sumably because of the large bulk conductivity. 



5. Field theory of diffusion 

In our discussion of WAL we have employed a semi- 
classical path picture and provided its connection to di- 
agrammatic calculations. Another theoretical approach 
that is commonly adapted in localization studies is the 
field theory approach. The field theory describing diffu- 
sion is the non- linear sigma model (NLctM)P21 The fields 
in the NLcrM live on a manifold that is determined by the 
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symmetry class. In some cases, these manifolds allow the 
presence of a topological term, that is a term that only 
depends on the topology of the field configurations, in 
the NLcrM. The localization properties depend strongly 
on the presence or absence of the topological terms (cf. 
section II B ) . 

To have a picture in mind, we provide a connection 
between the NLcrM and the semiclassical path approach. 
The propagators of the fields in the NLcrM are given by 
the diffusion and the cooperon, see figure [6] The non- 
linearity of the NLcrM corresponds to interactions be- 
tween the cooperons and diffusions. In particular, the 
four field interaction vertex corresponds exactly to the 
Hikami box, or the crossing point. Higher order inter- 
action vertices correspond to int erfer ence between paths 
with more than one self-crossing Ji2!2l 

The coupling constant in the NLcrM is related to the 
conductivity of the system. The WAL correction, as de- 
picted in figure [3j is the one-loop renormalization of this 
coupling constant. Higher order correction have been 
calculated and in the symplectic metal the next nonzero 
correction is only obtained at four- loop order (the calcu- 
lation has been carried out to five-loops) P^U This is one 
reason why the WAL correction ( 19 1 describes numerical 



simulations very well already at not too large conductiv- 
ity and system size (cf . section II B I. 



II. LONGITUDINAL CONDUCTIVITY OF 
TOPOLOGICAL INSULATOR SURFACES 

In graphene, the key observations confirming the Dirac 
nature of the charge carriers were the anomalous quan- 
tum Hall effect and the ambipolar field effect with the 
accompanying minimum conductivity!^^ Observation of 
these effects in TI's would constitute a convincing step to- 
wards verifying through transport experiments the Dirac 
dispersion of the surface states. Magnetic field and flux 
effects on quantum transport, including the quantum 
Hall effect, will be discussed in detail in the next two 
sections. In this section we focus on the zero field longi- 
tudinal conductivity and the absence of localization. 



A. Ambipolar field effect 

Observation of minimum conductivity and ambipolar 
field effect requires the conductance to be dominated by 
the surface. Progress towards that end has been made 
by reducing the bulk density by doping and using thin 
films. Several groups have reported minimum conductiv- 
ity and/or ambipolar field effect EHIhioZI As an example, 
figure [7] shows experimental data demonstrating both a 
minimum conductivity and an ambipolar Hall field ef- 
fect in Bi2Se3 thin films. As in graphene the minimum 
conductivity depends only weakly on temperature. The 
Hall density is linear in the gate voltage at low and high 
gate voltage, suggesting that the charge carriers are of 
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FIG. 7. Longitudinal resistivity of a Bi2Se3 thin film (about 
10 nm thick) obtained by mechanical exfoliation as a func- 
tion of back gate voltage at various temperatures (top). To 
reduce the bulk density the sample has been chemically p- 
doped with F4TCNQ molecules. The lower panel shows the 
carrier density as obtained by a Hall measurement. The lin- 
ear dependence of the Hall density at high and low gate volt- 
age indicates unipolar electron and hole transport. The inset 
shows an optical micrograph of the device with a scale bar of 
2 pm. Adapted from Ref. IT021 



only one type. Together these results indicate that the 
transport is through a surface state rather than being 
obtained from the bulk. 

In graphene the minimal conductivity is understood 
in terms of the physics of electron and hole puddles in- 
duced by disorder (for a review see Ref. IT0|) . The same 
is likely to hold true in TI transport, with nonlineari- 
ties in the spectrum away from the Dirac point playing a 
more important role than in graphene! 27 * 108 ! In fact, elec- 
tron hole p uddl es have already been observed in STM 



experiments!^' This is an important and current area 
of theoretical and experimental interest. It is, however, 
outside of the scope of this review which focuses more 
on the quantum correction to the conductivity. We will 
therefore not discuss it further, but rather briefly review 
the quantum theory of the longitudinal transport, single 
parameter scaling and the absence of localization. The 
results of these considerations will be of relevance to later 
discussion of WAL and quantum transport in WTI's, and 
is of fundamental interest by itself. 
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(a) massless Dirac model 



(b) random SO model 
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FIG. 8. Conductivity a at the Dirac point as a function of sys- 
tem size L obtained by a numerical solution of the Dirac equa- 
tion in presence of scalar disorder. The different data points 
are for different disorder strengths and correlation length. 
The system size has been rescaled by a mean free path to col- 
lapse the data on a single parameter scaling curve. At large 
conductivity and system size the scaling curves follows the 
logarithmic WAL dependence ( |19[ ) with the expected slope of 
1/n. The inset shows the scaling function /3 = dlog cr/dlog L 
demonstrating renormalization flow to the stable symplectic 
metal phase. (The factor of 4 in the conductivity scale is 
from the valley and spin degeneracy of graphene, and can be 
discarded for TPs). Adapted from Ref.lTTOl 



B. Absence of localization 

The surface of a strong TI is to pologically protected 
from Anderson localization.^HnSl Instead, disorder al- 
ways drives t he surface into the stable symplectic metal 
phaseP^- ° ! 112 !This even holds at the Dirac point, where in 
the clean case the density of states goes to zero and the 
condition for the perturbative calculation giving WAL 
(kp£ ^ 1) does not hold. The surface transport in TI's 
is therefore always diffusive at low enough temperature 
and characterized by WAL. 

Anderson localization^ generally happens in 2D elec- 
tronic systems when the dimensionless conductance g ~ 
1. In this regime, the perturbative calculations of sec- 
tion |I C| are not valid and one needs to rely on nu- 
merical simulations. In figure [8] we show numerical re- 
sults demonstrating the flow of the conductivity with 
increasin g syste m size L towards that of the symplec- 
tic metalP^- * 114 ! At large enough L the conductivity in- 
creases logarithmically with a slope 1/tt consistent with 
WAL p9|. The scaling function /3(a) = d\ncr/d\nL 



is strictly positive,! 110 * 114 ! as opposed to the normal 2D 
electron gas with s pin-or bit coupling which has a metal- 
insulator transition^ 13 ! 115 ! (and therefore a sign change in 
/9) at a ~ 1.4. 

The topological nature of the protection from localiza- 
tion can be understood by exploring the dependence of 




FIG. 9. The energy spectrum of a massless Dirac fermion 
(left) and a random spin-orbit coupled 2DEG (right) in pres- 
ence of disorder, as a function of boundary twist angle <p where 
ij}(x, y + W) = exp(i<p)ip(x, y). At ip = and ip = n the spec- 
trum is Kramers (doubly) degenerate. The Kramers pairs 
switch partners in the Dirac case but not in the 2DEG case. 
This is one manifestation of the Z2 topological classification 
in the symplectic class. Adapted from Ref. 11141 



the energy spectrum on the boundary condition! 111 ! 114 ! 
Specifically, figure [9] shows the sensitivity of the energy 
levels of a disordered Dirac fermion in a finite sample to a 
twist in the boundary condition ip(x, y+W) = e lip Tp{x, y). 
At ip = and tp — -k the system is time-reversal sym- 
metric and the spectrum is Kramers degenerate. Go- 
ing between these two special values the Kramers pairs 
exchange partners) 1 ^ This pair switch is a topological 
property and disorder can not change this behavior. The 
band width 8 V is always of the order of the mean level 
spacing A and by Thouless argument 1 - the conductance 
is at least g ~ <5 V /A <~ 1 and localization can not take 
place. Indeed, a localized state is unaffected by changes 
in the boundary conditions and its energy is independent 
of the boundary flux Lp and S v — > 0. In a conventional 
2DEG the Kramers partners do not switch pairs and lo- 
calization ensues for strong enough disorder. The switch 
versus no switch of Kramer's partners reflects the Z 2 na- 
ture of the topological classification of 3D time-reversal 
invariant TI'sPJ] 

At the field theory level, the difference between a single 
Dirac fermion and a 2DEG is understood to be due to the 
presence or absenc e of a topological term in the NLctM 
respectively! 112 ! 11 -^ The 3D TI is thus characterized by a 
topological term in the effective field theory of diffusion 
at the surface which is in one lower dimension. General- 
izing this notion to all dimensions and symmetry classes 
one can classify all possible topological insulators and su- 
perconductors by the possibility o f rea lizing a topological 
term in the corresponding NLcrMPSI 

In addition to the effects of generic disorder discussed 
here, several works, both theoretical and experim ental 
have explored the effects of local impurities jiiZHIHl 
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III. MAGNETOTRANSPORT IN CRYSTALS OF 
3D TOPOLOGICAL INSULATOR MATERIALS 

A. Shubnikov-de Haas oscillations and quantum 
Hall effect 

As a prelude to discussing the Aharonov-Bohm effect 
and other flux effects in the following section, we now dis- 
cuss standard magnetoconductance measurements on 3D 
topological insulators. Electrical transport in a 3D topo- 
logical insulator receives contributions from both bulk 
and surface states. A goal of much current research is 
to find 3D topological insulator materials that are truly 
insulating in the bulk, as defined for example by a resis- 
tivity that shows activated behavior at low temperature 
(i.e., diverges exponentially as exp(E g /ksT), where E g 
is the bandgap, up to power-law factors). This has been 
difficult to achieve, and in the few cases where a diver- 
gence of this form is seen, the effective bandgap is much 
smaller than the expected bandgap from either photoc- 
mission calculations or electric structure theoryP^Hl 

Magnetotransport measurements provide a natural 
means to distinguish between bulk and surface transport 
when the bulk retains a nonzero conductivity. Ideally, 
a sample can be made sufficiently clean that both bulk 
and surface contributions will show magnetoconductance 
oscillations, the Shubnikov-de Haas (SdH) oscillations. 
Such oscillations are periodic in inverse magnetic field 
with period 

A(1/B) = §J, (24) 

where Sf is the cross-sectional Fermi surface area trans- 
verse to the applied field. This is an area in fc-space, 
which has dimensions of inverse area, so the oscillation 
period is essentially an area divided by the flux quan- 
tum, up to a numerical factor. An ideal 3D topological 
insulator with one bulk Fermi surface and one surface 
Fermi surface will show two periods in the magnetocon- 
ductance oscillations, and these two oscillations will be 
distinguished by different and different angular de- 

pendences, assuming that the bulk Fermi surface is less 
two-dimensional than the surface one. 

In a strictly two-dimensional electron system, SdH os- 
cillations occur as a precursor of the quantum Hall effect. 
On a quantum Hall plateau the diagonal resistivity p xx 
drops to zero (note that the tensorial nature of p implies 
that the diagonal components of both resistivity and con- 
ductivity are zero). The locations of the Landau levels 
in a conventional 2DEG are at half-integer multiples of 
the cyclotron energy uj c oc B, E n — (n + l/2)huj c , with 
n = 0, 1, .... At fixed Fermi level Ep in the 2DEG, as is 
achieved if the 2DEG is in contact with a reservior, the 
condition for the nth Landau level to cross is E n = Ep. 
or 

(n + l/2)oc4- ( 25 ) 



So these features are equally spaced in 1/B, but with an 
offset of 1/2. In an STI this offset is absent because the 
spectrum is Dirac-like: including a Zeeman term for later 
reference, the Landau levels for a Dirac cone of velocity 
v occur at 

\E n \ = ^/(gpB/2) 2 + 2\nB\hev 2 , (26) 

where the first term is the Zeeman effect and the second 
gives an E ~ yB dependence. For zero Zeeman effect, 
the Landau level crossings occur at 

2nBhev 2 = E F 2 (27) 

so now n oc (1/B) with no offset, as promised. Hence a 
"Landau level index plot" ("fan diagram") can in principle 
distinguish between the case of a quadratic spectrum and 
a Dirac spectrum. 

In the SdH regime the same difference in offset shows 
up in the phase of the quantum oscillation. However, in 
SdH oscillations the diagonal resistivity and conductivity 
oscillate out of phase, while in the quantum Hall regime 
they move in-phase. For example, in the quantum Hall 
regime both are zero on a plateau as the resistivity and 
conductivity are purely off-diagonal tensors. A useful ref- 
erence for SdH oscillations is the book of Shoenbergp^ 
here we focus on the key differences between Dirac and 
ordinary (quadratic) fermions. A basic assumption in 
the following conventional theory is that the only effect 
of the magnetic field is that its coupling to the orbital 
degrees of freedom is essentially as a "probe": it does 
not modify the electronic structure except through this 
orbital effect. As high magnetic fields have to be used 
to observe oscillations, in part because the surface elec- 
tron density is high and mobility is not much higher than 
10 4 V/(cm 2 /s) in current samples even at low tempera- 
ture, this assumption needs to be considered carefully, 
especially in possible future materials where the surface 
topological state results from correlation phenomena that 
are sensitive to magnetic field. 

The two cases of SdH oscillations (normal 2DEG and 
Dirac fermions) can be obtained starting from the semi- 
classical quantization of energy levels: 

2ttpB 

S F (E) = —(n + 1 ) (28) 

where n is the Landau level index of the oscillation, Sp 
is the cross-sectional Fermi surface area at energy E and 
the shift 7 = 1/2 for a conventional 2DEG and 7 = 
for an STI. As emphasized by Mikitik and Sharlaij^HHH 
the shift is quantized to one value or the other depending 
only on the number of Dirac points; details of energetics 
are unimportant. While this semiclassical formula is only 
generally valid for large n, note that for 7 = 1/2 in a 
conventional 2DEG, we obtain for Landau level n 

S F (E)=nkp 2 ^ 2 ^ B (r l+ ^ 

=> E = hui c ( n+ - J (29) 
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consistent with the quantum Hall limit. For the Dirac 
case, we obtain at Landau level n 



2 2veB 
bp(-tj) = i^kp = — - — n 

=> E 2 = 2hnv 2 eB, 



(30) 



again consistent with the exact calculation. 

In practice the strength of the oscillations is reduced 
by both disorder and thermal fluctuations. The Lifshitz- 
Kosevich formula can be applied to estimate the strength 
of the oscillations incorporating this shift; in the simplest 
case (for additional corrections, see Ref. 11241 for experi- 
mental examples, see Ref. 11251 and references therein), 



AR tx RtRd cos 



~eB 



- 27T7 



(31) 



is the oscillatory component of the resistance R. Here 
Ro is the Dingle damping factor from disorder, 



Rd 



-tt/tu] c 



(32) 



where r is the scattering time, and Rt describes reduc- 
tion d ue to thermal smearing. One recent theoretical 
analysis^Sl of SdH experiments in Bi2Se3 and Bi2Tc3 
specifically for this phase 7 argues that all experiments 
are consistent with 7 = once effects of nonzero chemical 
potential and Zeeman energy are taken into account. Our 
discussion here has been in terms of flat infinite surfaces. 
The QHE effect on finite curved and multiple connected 
surfaces has been explored theoretically in Re fs. 1126111281 

There are several factors, some fundamental and some 
material-dependent, that complicate the observation of 
this physics in a 3D topological insulator. First, while 
the bulk of a topological insulator probably does serve 
to some extent as a reservoir of electrons for the surface 
state at a fixed chemical potential, this ignores surface 
charging effects and also ignores the dependence of bulk 
electron properties on applied magnetic field. Second, in 
transport the overall measured conductivity or resistivity 
always include a significant bulk contribution in current 
materials, although impressive progress has been made 
in reducing the bulk contribution. 

Another complication is that the g factor may be sig- 
nificantly larger than 2 (as large as 30-50); this is known 
to be the case for bulk electrons in Bi2Te3, Bi2Se3, and 
HgTe materials as a consequence of the same relativistic 
effects that lead to the TI behavior. Electronic struc- 
ture calculation of the g factor of the surface state is 
technically challenging, but STM experiments and most 
transport experiments are consistent with g < 10 so that 
the Zeeman effect is relatively weak compared to the or- 
bital effects. A "smoking gun" for the Zeeman effect is 
that the zeroth Landau level's energy is solely determined 
by the Zeeman effect, so observation of the zeroth level 
at high fields would allow a fairly direct measurement of 
the Zeeman effect, but this requires starting with a sam- 
ple whose chemical potential is close to the Dirac point. 



Some intrinsic interest of the Zeeman effect in the sur- 
face state is that it does not commute with the starting 
Hamiltonian, unlike in either graphene or a conventional 
2DEG (assuming spin-orbit coupling can be neglected). 
As a result, the Zeeman effect modifies the eigenstates, 
not just their energies, and opens up an energy gap even 
if the orbital field is excluded. 



B. Experimental observations of Shubnikov-de 
Haas oscillations and the quantum Hall effect 

Early observati ons of qu antum oscillations were domi- 
nated by the bulk! 129 " 131 ! From the SdH oscillations the 
anisotropic nature of the bulk Fermi surface is verified. 
This anisotropy also show s up in a ngular-dependent mag- 
netoresistance oscillations ! 130 * 13 ^ For these samples with 
large bulk conductance, a good characterization of these 
bulk contributions is required to extract any possible sur- 
face contribution. This may require very large fields 

Alternatively, in samples with larger sur- 
face co ntribution, cle ar 2D SdH oscillations are 
observed MMIMlmm The 2D nature of the SdH os- 
cillations is for example revealed by tilting the magnetic 
field. The 2D SdH only depends on the perpendicular 
component of the magnetic field, and therefore the posi- 
tion of resistance maxima varies as 1/ cos 9. An example 
of SdH oscillations in Bi2Te3 is shown in figure [10] 
The position of the maxima in the resistance can in 
principle be used to extract the value of the Berry phase 
through the shift 7 introduced in the last section. Most 
experiments have been interpreted to be consistent with 
a non-trivial Berry phase, though an accurate extraction 
of its value is complicated b y nonlinea rites of the spec- 
trum and a nonzero g factor.^ 139 * 140 ! An example of a 
Landau level fan diagram is given in figure |11| In many 
samples the doping level of the surface is time dependent 
due to aging effects. By using this time dependence it 
was possible to pr obe both electron and hole transport 
in the same samplef^ 

While there are several works that report SdH oscilla- 
tions, a fully developed QHE is rare th ough signatures of 
the QHE have been reported! 51 * 141 * 14 ^ In particular, fig- 
ure [12] shows the Hall conductance obtained in strained 
HgTe with plateaus at multiples of e 2 /h. The appearance 
of plateaus at even values of a xy / (e 2 / h) is interpreted to 
be a consequence of a different filling factor in the top 
and bottom surfaces. The longitudinal conductance does 
not go to zero at the plateaus, presumably due to the 
side surfaces remaining conductive. 

A quantized Hall effect and SdH oscillations with 
1 / cos 9 dependence have also been observed in highly 
doped samples, known to be dominated by the bulk,Q331 
and interpreted as evidence of layered bulk transport. 
As opposed to the surface data, these layered SdH 
oscillations where consistent with zero Berry phase. 
Finally we mention that related quantum oscillations 
have also been observed in other quantities, such as 
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FIG. 10. The resistivity derivative dp xx /dH as a function 
of the perpendicular component of the magnetic field 1/H± 
(upper panel) where H± = HcosO and 6 is the tilt angle 
between the magnetic field direction and the surface normal. 
The various curves correspond to different tilt angles. Non 
metallic samples of Bi2Te3 are obtained by growing with a 
weak compositional gradient. The lower panel shows the tilt 
angle dependence of the position of the third Landau level, 
which is consistent with the 1/ cos 6 expected from a 2D Fermi 
surface (solid line) . Adapted from Ref . | 
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IV. MAGNETIC FLUX EFFECTS ON 
QUANTUM TRANSPORT IN TOPOLOGICAL 
INSULATOR SURFACES 
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FIG. 11. Position of maxima and minmia of SdH os- 
cillations (Landau level fan diagram) in a 30 fim thick 
Bi1.5Sbo.5Te1.7Se1. 3 sample. The blue (left) data is obtained 
5 hours after cleaving when the surface chemical potential is 
such that the surface is p-type. After 720 hours the chemi- 
cal potential has shifted and the surface is now of n-type as 
demonstrated by the red (right) data. The solid lines are fit 
to theory of SdH taking into account nonlinearity of the spec- 
trum and is consistent with the expected Berry phase of it. 
Adapted from Ref. EH 
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FIG. 12. Hall conductivity as a function of magnetic field for 
a strained 70 mm thick HgTe sample, measured at 50 mK. 
The strain opens up a gap in the nominally semimetallic 3D 
HgTe, which is estimated theoretically to be of the order of 20 
meV. The inset shows the Hall and longitudinal resistances. 
While the Hall conductivity has some well developed plateus 
the longutidinal resistance does not go all the way to zero. 
Adapted from Ref. [5T] 



In this section we discuss two related effects of mag- 
netic flux on quantum transport, weak anti-localization 
and Aharonov-Bohm or Altshuler-Aronov-Spivak oscilla- 
tions. While these are often viewed as separate effects, 
they all arise from the AB phase of the wave function 
and quantum interference between paths. In particular, 
WAL and AAS oscillations are exactly the same phe- 
nomena except they are realized in samples with different 
topology (flat versus multiply connected) . When coupled 
with the Berry phase, AB oscillations of seemingly differ- 
ent nature can arise as discussed below. The distinction 



between the AB and the AAS oscillations is really that 
the former is realized when the particle motion is ballistic 
or pseudo-diffusive, while the latter are obtained in the 
diffusive state. An important difference between the two 
cases is that in the diffusive state, as we discussed in the 
introduction, one can not distinguish between a TI and a 
normal 2DEG. This is in principle possible with the AB 
oscillations. 
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FIG. 13. Magnetoconductance of a Bi2Se3 thin film grown 
on a Si(lll) substrate, for two different top gate voltages 
(left). The data has been fitted to the WAL expression ( |22[ ) 
to obtain the amplitude a and the phase coherence length 
Iq,. (The definition of a in this plot is of opposite sign to the 
one adopted in this review). The plot on the right shows the 
extracted values of a and as a function of the top gate 
voltage. A large negative gate voltage induces a depletion 
layer between the top surface and the bulk channel which 
then independently contribute to the WAL. Adapted from 
Ref.[T53l 



A. Weak anti-localization in thin films 

Since the surface state is always driven into the sym- 
plectic metal phase, one expects to observe weak anti- 
localization. WAL can both be observed as a negative 
magnetoconductance ( |22| and as a logarithmic tempera- 
ture enhancement of the conductivity (19 1. We will dis- 
cuss the two approaches separately in the following. To 
avoid having the WAL signal being completely masked 
by bulk effects most experiments t o date wo rk with thin 
films. These are mostly Bi 2 Se 3 film d 95 ' 96 l 147 ' M3but stud- 
ies with Bi 2 Te^Hland Bi^Se^Tei.^^have also been 
reported. To further reduce the bulk conductivity the 
films are som etimes add itionally doped, for example with 
CaJ^S PqUSU or c u li55| or |- ne p erm j level is moved into 
the bulk gap with the help of a gate voltage on a back 
gate^EHEa or a top gatePH The effect of magnetic 
doping and the resulting crossover to weak localization 
(WL) has also been studied, both experimentalljff^ and 
theoretically.^ 



Magnetic field dependence of weak anti-localization 



Most experiments observe a negative magnetoconduc- 



tance consistent with ( 22 1 but with a p refactor that varies 
roughly in the range 0.3 < a < i.2j 95 | 96 | i47ffl gBIIggIIBa 
The value of a can be tuned by as much as a factor of 
two in a single sample by varying gate voltage^- 48 * 1 5 ° l 1 53 l 
(see figure 13 1. Both the broad range of values obtained 
for a and its gate voltage tunability can be understood 
by carefully taking into account the bulk contribution 
and the coupling between surface and bulk. Some of this 
physics was discussed theoretically in Ref . [75] and later 



in detail in Ref. 1751 We give a simplified discussion of 
the results of Ref. [TS] here and refer to the original works 
for further details. 

When the thickness of the film t <C £</>, the bulk is ef- 
fectively two dimensional and contributes 2D quantum 
correction of its own. Generally, due to the strong spin- 
orbit coupling in TI materials, this is also a WAL of the 
form (22) with a = 1/2. However, if the bulk is only 



weakly doped and intervalley and spin scattering in the 
bulk is negligible, a WL with a = — 1 is obtained instead. 
Roughly speaking, this is because the spin of the eigen- 
states close to the band edge depends only weakly on the 
momentum. The Berry phase from the spin rotation is 
thus absent. It is likely that in most current samples the 
bulk is in a parameter range where WAL with a = 1/2 
is expected. Strong coupling with the surface state will 
also generally give rise to WAL, even if one expects WL 
from the isolated bulk. In the following we assume the 
bulk is in the WAL regime. We note though that WL has 
been observed in ultrathin films (4-5 quintuple layers) 
in Bi 2 Se 3 nlmsPfil and (Bi :c Sb 1 _ a; ) 2 Te 3 films. 1 ™ It does 
not seem likely that these observations are explained by 
the bulk contribution. In particular, in the former study 
the presence of WL is substrate dependent, suggesting 
that details of the hybridization of the top and bottom 
surface and their environment plays an important role. 

If the surface and bulk are completely decoupled their 
contribution to WAL will simply add and a magnetocon- 
ductance (22 1 with a = 1 is expected. In contrast, if the 



bulk and surface are strongly coupled, r s b <C with 
the surface to bulk tunneling rate, they act as a single 
channel with a = 1/2. With intermediate surface-bulk 
coupling a value of a somewhere between the two limiting 
values is expected. Top and bottom surface can also con- 
tribute independently and have different coupling with 
the bulk and different phase coherence lengths. This ex- 
plains the observed range of values for a obtained in the 
experiments. The gate voltage dependence of a is likely a 
consequence of a variation of the surface-bulk coupling. 
The gate voltage induces a depletion region below the 
surface, reduci ng the tunnel coupling of the surface to 
the bulkPmSl 

Our discussion in this section (and most analysis of cur- 
rent experiments) has been simplistic in the sense that 
it assumes that the observed WAL correction to the con- 
ductance can still be described by (22 1, just with a renor- 
malized value of a. While there are limits in which this 
is justified,^ the bulk, top and bottom surface generally 
have different phase coherence lengths and the WAL is 
described rather by a sum of terms of the form ( 22 1 . Nev- 



ertheless, the good qualitative agreement between theory 
and experiment can be taken as a convincing evidence for 
surface transport. However, since WAL is a characteris- 
tic of the symplectic metal that is independent of the 
topological properties of the system realizing the phase, 
a different type of experiments is needed to probe the 
topological properties of TI surface states through trans- 
port. 



1G 



2. Temperature dependence of weak anti-localization and 
the effect of electron- electron interactions 

In contrast to the good agreement between magne- 
toconductance experiments and the theory of the sym- 
plectic metal, the temperature dependence of the con- 
ductance does not show WAL following (19). In- 
stead a WL behav ior, that is (19 1 with a negative, is 
observedP^ ti 151 * 1 -^ This may possibly be explained by 
taking electron-electron interactions into account, which 
we have so far ignored in our discussion. Electron- 
electron interactions also give a logarithmic correction 
to the con duct ance, the so called Altshuler-Aronov (AA) 
correctionP^ This correction is relatively independent of 
a magnetic field and can therefore in principle be differ- 
entiated from the WAL correction by the application of 
a magnetic field, as has been observed experimentally. 150 

These experimental observations are a least qualita- 
tively consistent with the AA scenario, but convincing 
quantitative comparison is lacking. Other explanations 
for these observations ha ve not been ruled out, or even 
explored theoretically. 162 If, however, these are really sig- 
natures of electron-electron interactions, it opens up the 
possibility of observing interesting correlation physics. 
For example, it has been suggested that the combination 
of the localization effect of interactions and the topolog- 
ical protection of the surface will drive the surface into a 
stable critic al st ate, characterized by a universal value of 
conductivity.^ 



B. Aharonov-Bohm effect in topological insulator 
nanowires 



An ideal TI nanowire, with an insulating bulk, can be 
thought of as a hollow metallic cylinder. As explained in 
the introduction, a magnetic flux along the length of the 
wire leads to periodic oscillations in the conductance as a 
function of flux with period h/2e: the AAS oscillations. 
These oscillations are simply another manifestation of 
WAL as they are a consequence of quantum interference 
between clockwise and anti-clockwise circulating paths. 
As such, the AAS oscillations are characteristic of the 
symplectic metal and can not directly probe the topo- 
logical properties of the surface state. Observation of 
these oscillations is however a strong indication of sur- 
face transport. 

The AAS oscillations are expected when the transport 
is diffusive. Oscillations with period h/e, twice that of 
the AAS oscillations, can be realized when transport is 
close to being ballistic, or when the Fermi level is at or 
close to the Dirac point. The basic physics behind these 
latter oscillations, to be discussed in the next section, 
is the Berry phase and the perfectly transmitted mode. 
Observation of these oscillation would constitute an in- 
direct measurement of the surface state Dirac fermion's 
Berry phase. 
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FIG. 14. A schematic drawing of the energy spectrum of a 
TI nanowire for two different values of flux along the wire. 
The horizontal axis k is the momentum along the length of 
the wire. The different branches of the spectrum correspond 
to different values of the transverse momentum q n = 2irn/W 
with W the circumference. All the lines are therefore doubly 
degenerate, corresponding to clockwise and anti-clockwise cir- 
culation around the wire, expect for the linear curves which 
correspond to q — and are nondegenerate. Consequently 
the number of transmission modes at flux cf> — (j>o/2 is odd 
but even ad (j> — 0. Adapted from Ref. 11661 



1. Berry phase and the perfectly transmitted mode 

A TI nanowire differs from a conventional quantum 
wire in that it can host an odd number of transmission 
modes and therefore a perfectly transmitted mode with 
conductance of e 2 /h. Due to the Berry phase the per- 
fectly transmitted mode is realized at a flux thro ugh the 
wire equal to a half integral number of flux quantaJ^allSSl 
If all other modes in the wire give a negligible contri- 
bution to the conductance the conductance will oscillate 
with a period of h/e with a minimum at flux <fi = 0.l 166 l 167 l 

To demonstrate this in an explicit model we consider 
the effective theory of the surface state (some of the re- 
sults belo w have also been obtained with a 3D lattice 
model 167 ), described by a single Dirac cone 



H 



vp ■ a. 



(33) 



In addition to the kinetic term a curvature induced spin 
connection term, that describes how the spin rota tes as 
it moves along the surface, is generally present P^H For a 
cylindrical surface this term can be completely absorbed 
into the boundary condition 



i/)(x,y + W) = e i7r ip(x,y), 



(34) 



which is now anti-periodic. Here x is the coordinate along 
the wire and parallel to the flux, y is the transverse coor- 
dinate and W is the circumference. The anti-periodicity 
is a consequence of the fact that the spin lies in the tan- 
gent plane to the surface and therefore rotates by an angle 
of 2tt when looping around the surface. Because of the 
antiperiodic boundary condition the spectrum is gapped 
as schematically shown in figure [14] Within the same 
approximation, the effect of applying a flux along the 
wire is simply to change the boundary condition by the 
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FIG. 15. Two terminal conductance of a TI nanowire as 
a function of Fermi energy e = hvkF, for various disorder 
strengths Kq. For each value of Kq the conductance is plot- 
ted for three different values of flux <f). When the conductance 
at <f> = and (j> — 4>o/2 is the same the magnetoconductance 
oscillates with a period of h/2e (AAS oscillations). This hap- 
pens at moderate to large disorder strength and at doping 
away from the Dirac point. Magnetoscillations with period 
h/e are otherwise obtained. At the Dirac point the conduc- 
tance always has a maximum at flux (f> — 4>o/2 but away from 
the Dirac point a maximum at either <j> = or cj> = <j>o/2 is 
obtained depending on the doping level. The data is obtained 
by numerically solving the Dirac equation in presence of dis- 
order for a wire with W = 100£ and L — 200£ where £ is the 
correlation length of the disorder. The oscillations as a func- 
tion of Fermi energy away from the Dirac point at Kq = 0.2 
have a period which corresponds to the mean level spacing 
A = Hv/W. Adapted from Ref.lTrJol 



Aharanov-Bohm phase 

4>(x, y + W) = e*( 2 ^o+?0 ^ y y (35) 



A flux of half a flux quantum, <j> = 4>o/2, cancels the Berry 
phase and the spectrum becomes gapless (figure 14 1. 
Crucially, the number of modes is now odd and there is 
necessarily at least one perfectly transmitted mode with 



conductance e /h. 

Figure [15] shows the conductance as a function of dop- 
ing as obtained by a numerical simulation)!^ The simula- 
tion includes a Gaussian disorder potential V(r) that sat- 
isfies (V(r)V(r')) = K (hv) 2 exp[(r - r')7(2£ 2 )]/27r£ 2 . 
Ko is a dimensionless measure of the disorder strength 
and £ is the length scale characterizing the disorder (typ- 
ical size of electron-hole puddles). For each value of Kq 
the conductance is given at three values of the flux, which 
determine the period of possible oscillations. Three dif- 
ferent regimes are observed: i) If disorder is large enough 
such that transport is diffusive the AAS oscillations dom- 
inate. The conductance is the same at <f) = and 
4> = 0o/2- ii) At the Dirac point conductance is dom- 
inated by the perfectly transmitted mode. The conduc- 
tance oscillates with a period of h/e and has a minimum 



FIG. 16. Aharonov-Bohm oscillations in the resistance as a 
function of magnetic field of a Bi2Se3 nanoribbon of a rect- 
angular cross section (120 nm x 50 nm) and length ~ 2 yum. 
The left inset shows the position of resistance minima and 
the right inset the Fourier transform of the derivative dR/dB, 
with a strong peak at a magnetic field corresponding to a flux 
of h/e and a weaker peak at flux h/2e. The resistance has a 
minimum at zero flux and the relatively sharp increase of re- 
sistance around zero flux is indicative of WAL, possibly from 
the bulk which density is nonzero. Adapted from Ref. 168 



at zero flux, iii) At weak disorder and large enough dop- 
ing, such that the conductance is no longer dominated 
by the perfectly transmitted mode, the conductance os- 
cillates with a period of h/e. In this regime the conduc- 
tance can have either a maximum or a minimum at zero 
flux, depending on the doping level. The conductance 
oscillates as a function of doping with a period equal to 
the level spacing A = hv2n/W . These oscillations are 
smeared by temperature when ksT sa A and the con- 
ductance at (j) = and <j) = 0o is equal. The conductance 
then oscillates with a period of h/2e. 



Experimental observations of Aharonov-Bohm 
oscillations 



Aharonov-Bohm oscillations have been experimentally 
observ ed in n anowires (ribbons, plates) of Bi2Se3p^ 
Bi 2 Te 3 p21iZ0] an d /3-Ag 2 TeP Related periodic oscil- 
lations h ave also been observed in arrays of ringlike 
structures J^D In the wires, the magnetoconductance 
shows periodic oscillations superimposed on the UCF (see 
figure [16] for example data). Fourier transform reveals a 
strong h/e period and a weaker h/2e period. 

All the samples are doped away from the Dirac point 
which in Bi2Te3 is buried in the valence band. In Ref. 11691 
a back gate was used to tune the Fermi level and reduce 
the bulk contribution to the conductance. The AB os- 
cillations become more pronounced as the surface con- 
tribution is enhanced. If these observations are to be 
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consistent with the theory of the last section, the sam- 
ples would need to be in the weak disorder limit. As a 
necessary condition the mean level spacing must be large 
compared to the temperature broadening. From the cir- 
cumference of the samples one estimates A ~ 10A" — 70K 
while the measurement temperature is usually T ~ 2K. 
In Ref. 11701 a crossover from h/e to h/2e oscillations is 
observed at T « 5A for a sample with A m 10A. 

In the weak disorder regime the conductance at zero 
flux will be either maximum or minimum depending on 
the sample. In Refs. 155111681 and !169l the conductance has 
a maximum while in Ref. 170 it has a minimum. A change 
in the sign of the magnetoconductance with varying gate 
voltage was not observed}^' though a systematic explo- 
ration of this effect was not undertaken. A clear "smoking 
gun" signature of Dirac fermion transport in this regime 
would be the observation of the gate voltage oscillations 
accompanied by a ir shift in the magnetoconductance. 

Taken together these results are consistent with the 
theory of doped and weakly disordered TI wires. A more 
systematic exploration of doping and disorder strength 
dependence of the AB oscillations is needed to firmly 
confirm this picture. In particular, the robust h/e oscil- 
lations at the Dirac point and the accompanying indirect 
measurement of the Berry phase is still to be observed. 



V. RELATED PROBLEMS AND FUTURE 
DIRECTIONS 

We end this review by briefly exploring a couple of 
related topics. These arc quantum transport in WTI and 
the Josephson effect in TI surfaces. 



Transport at the surface of a weak topological 
insulator 



So far we have mainly focused on the transport prop- 
erties of strong topological insulators (denoted in this 
section with STI to clearly differentiate them from the 
WTI). STI's have an odd number of Dirac cones at their 
surface, while weak TI's have an even number of Dirac 
cones. A WTI may be realized e.g. in layered semicon- 
ductors, such as KHgSbJ^l Topological crystalline insu- 
lators which surface electronic structure is also described 
by an e ven num ber of Dirac cones, have recently been 
proposed 173 * 174 ! and observed Ji^iZZI In addition a thin 
film of STI can be considered to be a WTI. 

A WTI can be thought of as being obtained by stacking 
2D QSHE layers. The helical edges states of the QSHE 
couple together to form the surface state of the WTI. By 
construction this surface state will not be found on every 
face of the sample. If the number of layers is odd the 
number of propagating modes in the surface is also odd 
and the surface hosts a perfectly transmitted mode. The 
conductance is at least e 2 /h and the surface can not local- 
ize. The energy levels have a similar switching behavior 



as m 

STlPS Instead of localization the surface of a WTI 
flows, under rather general conditions to be discussed be- 
low, into the symplectic metalf^ l n this sense the WTI 
behaves very much like an STI. The effect of interactions 
on WTI surface states was studied in Ref. 11801 

The surface state of a WTI can be modeled by two 
Dirac cones 



H = vt ® a ■ p + V(r). 



(36) 



Here p is the momentum operator and the Pauli matrices 
a and r act in spin and valley space respectively. The 
disorder potential V can be written as 



aB 



(37) 



1, t z £g> 1 and 



Of all the possible terms, only 1 ® 1 , - 
T y ® cr 7 with 7 = x, y, z do not break the time-reversal 
symmetry which is given by T = 1 ® ia y C with C the 
complex conjugate, and satisfies T 2 = — 1. Of these, 
only T y <g> <J Z opens up a gap in the energy spectrum and 
we therefore refer to the average of this term as the mass 
m = (r y <£> o z ). The disorder terms are independently 
distributed with (5V a p(Y)5V a p{v')) = gK(r — r') where 
J drK(r) = 1. 

For a random potential the mass is zero. A nonzero 
mass requires the surface potential to be commensurate 
with an even number of unit cells, and is therefore nec- 
essarily zero for an odd number of layers. In this way 
the mass of the 2D theory connects with the even-odd 
argument above. Indeed, in the absence of mass dis- 
order drives the surface state into t he s ymplectic metal 
following a single parameter scalingP^ Figure 17 shows 
the result of a numerical simulation of ( 36 1 for various 



disorder strengths and doping levels including the Dirac 
point. The system size, horizontal axis, has been scaled 
by the mean free path to reveal a collapse of the raw 
data (inset) onto the scaling curve. At large system sizes 
the conductivity follows the logarithmic dependence of 
the WAL (19 1 with a slope of 1/V, just as in an STI. 
This is expected since the two Dirac cones are strongly 
mixed by the disorder and was verified with an explicit 
diagrammatic calculation in Ref. 11781 Note that while 
the WTI is similar to graphene in the sense that both 
systems have two Dirac cones, they differ in the sign 
of the square of the time-reversal operator. For a WTI 
T 2 = — 1 but for graphene T = 1. Intervalley coupling 
in graphene therefore lead s to weak localization^ and 
eventually localization. 181 

If m is nonzero, the surface can be localized by disor- 
der. A metal insulator transition is obtained at a con- 
ductiv ity a s s 1.4, just as in a regular spin-orbit coupled 
2DEGP^3The conductance distribution at the transi- 
tion is also the same as in the 2DEGP^This is consistent 
with the two systems being described by the same effec- 
tive field theory of diffusion, t he s ymplectic NLctM in the 
absence of a topological termini The fact that localiza- 
tion is avoided at m = suggest that the full scaling 
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FIG. 17. Conductance as a function of scaled system size 
for various chemical potentials n and disorder strengths g at 
zero mass m = 0. The inset shows the raw data before scal- 
ing. The data is obtained by a numerical simulation of the 
Dirac equation ( |36| |. At large system size the conductivity 
approaches the WAL form (19 1 with slope l/n (grey solid 
line). The conductivity follows single parameter scaling with 
a strictly positive beta function with an absence of Ander- 
son localization, just as in the strong topological insulator. 
Adapted from Ref. 11791 



flow may be described by a two parameter scaling, as in 
the quantum Hall effectP^ 3 -'. This possibility has been ex- 
plored in Ref. 11791 If the disorder becomes strong enough 
that nonlinear in momentum terms need to be added to 



the Dirac theory (36 1 localization is expected. 



The NLerM of a WTI does not have the topological 
term that protects the STI from localization.^^ How- 
ever, it has recently been argued that for localization to 
happen, a proliferation of vortex like topological defects 
in the field configurations is necessaryP^ Such vortex 
proliferation is forbidden by average discrete symmetries 
which are equivalent to having a zero mass, consistent 
with the observed absence of localization. 



a flux quantum is of course the magnetic flux associated 
with an elementary vortex in a charge-2e superconductor, 
so the setup is quite natural. 

The answer is that the metallic mode becomes gapped 
along the length of the vortex, but there are zero -energy 
Majorana excitations at the two ends of the vortexJ22Hi£8| 
(i.e., where it leaves the topological insulator). Their 
emergence in this system can be understood in two ways. 
Looking at the surface state of the topological insulator, 
the proximity effect converts the single-sheet Fermi sur- 
face into a "topological superconductor" wit h one M a- 
jorana fermion excitation in each vortex 

corePHHni An 

ordinary spinless fermion excitation can be broken up 
into two Majorana fermions, which are their own Hermi- 
tian conjugates. This topological superconducting layer 
is like a time-reversal symmetric version of the 2D "p+ip" 
superconduc tor known to support Majorana fermions in 
vortex cores An alternate picture of where the Majo- 
rana fermi ons c ome from is as tail states of the gapped 
vortex core. 192 For recent reviews on realizations of Ma- 



jorana modes in condensed matter systems, see Refs. 
and HH 

As a first step in the direction towards the observa- 
tion of the Majorana mode, and of intrinsic interest by 
itself, proximity effect and a Josep hson cu rrent through 
TI surfaces has been demonstratedP^M200| B t n Shapiro 
steps in the presence of microwave irradiation^^ a nd 
Fraunhofer interference pattern induced by a magnetic 
fi e ij i96 | i 98| fr ave been observed. The temperature depen- 
dence of the data in Ref. 11971 was more consistent with 
ballistic transport than diffusive transport, and the su- 
percurrent was found to be mainly carried by the surface 
state, even if the normal state transport was heavily in- 
fluenced by the bulk states. This servers as an additional 
motivation for Josephson physics as a probe of the unique 
properties of the normal state physics which is the focus 
of this article. In addition , no t all aspects of the observed 
data are fully understoocH^SI which should spur on fur- 
ther theoretical and experimental efforts in this line of 
work. 



B. Superconductivity and topological insulators 

Another exciting direction that merits a review article 
of its own is how the proximity effect from a conventional 
s-wave superconductor has unusual effects on topologi- 
cal insulators and related electronic structures such as 
strongly spin-orbit coupled semiconductors in a magnetic 
field. Here we mention only some aspects of this problem 
that relate to the discussion above of the Aharonov-Bohm 
effect. In a non-superconducting nanowire made from a 
topological insulator, there is a protected mo de when one 
half of a flux quantum pierces the nanowire! 163 * 165 ! The 
same metallic mode is there in the inverse problem where 
space is filled by a 3D topological insulator that has a 
cylindrical hole in it, pierced by one half of a flux quan- 
tum. What happens to this system when the 3D topo- 
logical insulator is driven superconducting? One half of 



VI. SUMMARY 

In this report we have reviewed recent theoretical and 
experimental progress towards the verification of Dirac 
fermion surface states in topological insulators, through 
various magnetotransport effects. We have adapted a 
positive outlook on the status of the field. While there is 
certainly a significant number of transport experiments 
that show no reliable signs of surface transport for vari- 
ous reasons, the evidence for Dirac fermion transport and 
the accompanying Berry phase is mounting. Material im- 
provements are rapid and the first steps in combining TI's 
with superconductors are being taken. TI transport ther- 
fore has the promise to be fertile grounds for exploring 
and realizing interesting physics. 

At the theoretical level, many of the transport proper- 
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ties of TI surfaces result from disorder driving it into the 
symplectic metal phase. The symplectic metal is charac- 
terized for example by weak anti-localization. However, 
since this phase is also realized in traditional two dimen- 
sional electron gas in the presence of spin-orbit coupling, 
the nontrivial Berry phase of the Dirac fermion is not 
easily realized in this limit. The Berry phase is revealed 
in principle in quantum oscillations such as Shubnikov-de 
Haas oscillations and, in TI nanowires, Aharonov-Bohm 
oscillations. We have discussed these effects in detail in 
the main text. 

Signatures of surface transport in large bulk crystals 
are weak if not absent. This is due to significant bulk con- 
duction, coming from unavoidable bulk doping. Instead, 
the most convincing data is obtained in thin films, where 
the surface can have a significant contribution to the to- 
tal conductivity. Observation of weak anti-localization 
fitting that expected from the symplectic metal, and 
2D SdH oscillations strongly suggest surface transport, 
though one should keep in mind that the bulk of thin 
films is often effectively 2D also. Comparison between 
theory, taking into account the bulk, and experiments 
suggests the presence of at least two transport channels, 
bulk and surface. Aharonov-Bohm oscillations have been 
observed in TI nanowires. These oscillations can in prin- 
ciple reveal the presence of the Berry phase, but this 



has not been attained with current data. SdH oscilla- 
tions are in most cases consistent with a nontrivial Berry 
phase, though the data can be difficult to interpret due to 
various correction. Finally, the experimentally observed 
temperature dependence of WAL is seemingly inconsis- 
tent with the symplectic metal, possibly pointing towards 
important effects of electron-electron interactions. 

In this review we have focused on the surface theory 
and mostly ignored the bulk. This reflects our view that 
eventually insulating bulks will be realized, and that the 
surface transport contains the more fundamental physics. 
Lastly, we mention that in addition to the reviews already 
cited in the introducti on! 1 *"^ a couple of focused reviews 
have appeared recently ^ 01 * 202 ! which discuss transport of 
TPs from a different perspective. 
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